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Abstract. We present an algorithm which computes a planar 2-spanner from an Unit Disk Graph when 
the node density is sufficient. The communication complexity in terms of number of node's identifier 
sent by the algorithm is 6n, while the computational complexity is 0(nA), with A the maximum 
degree of the communication graph. Furthermore, we present a simple and efficient routing algorithm 
dedicated to the computed graph. 

Last but not least, using traditional Euclidean coordinates, our algorithm needs the broadcast of as few 
as 3n node's identifiers. Under the hypothesis of sufficient node density, no broadcast at all is needed, 

iq reducing the previous best known complexity of an algorithm to compute a planar spanner of an Unit 

rsj Disk Graph which was of 5n broadcasts. 

1 Introduction 

In many problems on networks, among which the problem of message routing [1-3], it is useful to know a 
i— i planar subgraph of the communication graph. Although, all planar subgraphs are not equally interesting: 
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an usual requirement is that the length of a path between two nodes is not too much longer in the planar 
subgraph than in the original graph. 

In this paper, we propose a distributed and simple way to compute such a planar subgraph of a unit disk 
l/~j graph when the nodes of the communication graph are localized using the virtual raw anchors coordinate 

t-H system [4], instead of the stronger hypothesis of having the nodes localized in a classical 2D coordinate 

^f~\ system. 

We further propose a simple, efficient and light routing algorithm that is dedicated to the constructed 
graph. This last contribution is related to the conjecture (partially solved in [5] and fully solved [6, 7]) that 
a 3-connected graph accepts an embedding such that the greedy routing algorithm 1 is guaranteed to deliver 
any message. Indeed, we propose a routing algorithm {similar to the greedy routing), which ensures message 
delivery in the unit disk graph induced by any set of nodes, under some connectivity assumptions. 

X 

1.1 Related work 

Planar graph and poset dimension [8] it is proved that a graph G = (V, E) is planar if and only if 
it has order-dimension at most three, where a graph has dimension d if and only if there exists a sequence 
<i, . . . , <d of total orders on V whose intersection is empty, and such that for each edge (xy) € E and for 
each z G V \ {x, y}, there is at least one order <j in the sequence such that x <j z and y <j z. 

It means that any three total orders whose intersection is empty induce a planar graph, this graph being 
the subgraph of the complete graph obtained by keeping only the edges that satisfy the second condition. 
We will refer to this graph as the Schnyder's graph of the three total order, and we note it G< "^ < 3 or G s 
for short. 



1 Given a distance among the nodes, for instance the Euclidean distance, we call the greedy routing algorithm, the 
algorithm which consists for a source x to forward a message to the node that is the closest to the destination y in 
that distance 



Planar spanner much work have been dedicated to the construction of planar subgraphs. One of the 
first planarization technique is the use of Gabriel graphs, but, if connectivity is preserved, an edge may be 
replaced by a path of unbounded length [9], whereas, as we mentioned previously, we want to avoid this. 
The following well known definition catches the type of subgraphs we are interested in: given a graph G and 
a subgraph G' , we say that G' is a k-spanner if, for all pairs of nodes x, y, a shortest path from a node x 
to another node y in G' , is at most k times longer than a shortest path between these two nodes in G; the 
factor k is called the stretch. If G is the complete graph, we further say that G' is a geometric fc-spanner. 

The Delaunay's triangulation of a set of vertices V is a planar geometric spanner. Its stretch factor is 
upper bounded by 1.998 [10], and lower bounded strictly by 7r/2 [11], the exact stretch being unknown. In 
[12], the authors efficiently construct a planar 2.5-spanner of UDG that contains the edges of length 1 of 
the Delaunay's triangulation of a set of nodes V. The complexity of this construction is improved in [13], in 
which an algorithm needing 5 broadcast is proposed. An other construction of spanner of Unit Disk Graph 
is proposed in [14] with stretch > 2. Other spanners exist, in particular a way to construct a 2-spanncr from 
a complete graph is proposed in [15]. Interestingly, it is shown in [16] that three different constructions lead 
to the same planar geometric 2-spanner. These three constructions are the half-#6-g ra phs, the triangular- 
distance Delaunay triangulation (TD Delaunay graphs) and the geodesic embeddings. In [17], the authors 
further propose a planar spanner with bounded degree. We refer the interested reader to the recent survey 
of Bose and Smid [18]. 

Greedy embedding When one consider building a spanner, one usually does not focus on preserving easy 
routing properties such as greedy routing. It means that even if the greedy routing algorithm delivers a 
message to the destination in the original communication graph, it has no guarantee to succeed in an usual 
spanner. Preserving such a property would be of great interest. 

This problem is related to the following conjecture [19]: given a 3-connected graph, does an embedding 
exist such that the greedy routing algorithm is guaranteed to deliver any message ? 

In [5], it is proved that the conjecture is true if the graph is a plane triangulation by using Schnyder's 
characterization of planar graphs [8]. Later, the conjecture was proved for every 3-connected graphs in [6, 7]. 

1.2 Summary of results 

VRAC For all the results presented in this paper, we assume that the nodes are localized using the virtual 
raw anchors coordinate system (VRAC, [4]), or a simple variant of it. Supposing that the nodes are localized 
in these coordinate systems, is a strictly weaker hypothesis than the hypothesis that the nodes are localized in 
a traditional 2D coordinate system. Indeed, if the nodes are localized in a traditional 2D coordinate system, 
it is possible to compute their coordinates in the virtual raw anchors coordinate systems, while the converse 
is impossible. Furthermore, this coordinate system is expected to be easier to implement in practice. 

Planar subgraph and spanner As mentioned previously, a graph of order-dimension 3 is planar [8]. Hence, 
to planarize a graph, it is sufficient to select edges that correspond to three total orders. This technique has 
already been used, and, for instance, the half-0g graph mentioned in the previous section can be constructed 
along this line. However, several issues appears. First, this technique may need important computations. 
Second, the three orders being total, the computation may not be feasible locally. Third, the computed 
planar graph is a subgraph of the complete graph, and may not be a subgraph of the communication graph. 
In this paper, we address these issues when the communication graph is a Unit Disk Graph. To do so, we 
propose three total orders based on the VRAC coordinates using which we can construct a planar graph 
such that if the node density is high enough, 

1. it needs only comparison (no other operations of any type), 

2. it is a 2-spanncr 

3. constructing it requires to broadcast at most 6n nodes identifiers. 



In particular, our result improves the result of [12] by construction a spanner with stretch factor 2 versus 
2.5. Plus, using the VRAC coordinates, our algorithm induces the broadcasts of at most 6n node's identifiers 
(excluding the one needed for the neighborhood discovery), and has computational complexity O(nA), for 
A the maximum degree of G. Furthermore, using traditional Euclidean coordinates, it needs the broadcast 
of as few as 3n node's identifiers that can be done in a single communication round. As when the density is 
high enough, the constructed graph is a planar 2-spanner, our work answers the open problem number 22 of 
[18J. 

If this work is inspired by the paper of Schnydcr [8] , we stress out that the constructed graphs are not 
necessarily subgraphs of G , the planar graph induced by the three total orders as when following Schnyder's 
theory. 

In more details, in Section 3 we construct a first subgraph G from G by using only each node's neighbors. 
If we only keep the edges of length at most 2r/-\/5 « 0.8944r, the subgraph G is planar. When the node 
density is too small, the obtained graph may not be connected, and, to avoid this, we introduce virtual edges 
(edges that are not edges of the connexion graph). In Section 4, changing slightly the VRAC coordinates 
system, we prove that Gi is a subgraph of the half-6> 6 graph which is equals to G s . It implies that Gi is planar, 
but it gives no result on its stretch. Nonetheless, we prove that it verifies 1) the length of a shortest path in Gi 
is at most twice the length of a shortest path in G, and 2) a virtual edge corresponds to a path of two edges in 
G. When the node density increases, the virtual edges disappear, hence, Gi is a planar 2-spanner of G when 
the node density is high enough. All these results hold even when the constructed graph is not equal to the 
half-06 graph. Finally, using the VRAC coordinates, our algorithm needs two round of communications and 
induces the broadcasts of 6n node's identifiers on top of the one needed for the neighborhood discovery. Using 
Euclidean coordinates, we can reduce this to the broadcasts 3n node's identifiers that can be performed in 
a single round of communication, and no messages are exchanged at all when the density is high enough. 



Routing In [5], it is proved that a plane triangulation has an embedding in which the greedy routing 
algorithm is guaranteed to deliver any message. In our work, we assume given the embedding in the plane, 
so, it means that, instead of choosing an embedding for the nodes, we look at the dual problem, that 
is designing a routing algorithm (as close as possible to the greedy routing algorithm) which guarantees 
delivery. 



2 The model 

2.1 Communication model 

We consider a wireless network in which two nodes can communicate if they are at distance at most r, the 
communication radius. We can normalize the distances so that r — 1, in which case we have Unit Disk 
Graphs (UDG). However, we will keep mentioning r, as we believe it carries useful information. The use of 
the UDG model for the communication links is subject to caution from a practical point of view. We quickly 
mention the recent paper [20] that discusses how protocols that are proved valid under the UDG model can 
be turned to valid protocols in the more realistic SINR model. Another way of extending the results of this 
paper to more general communication models is to use basic properties of such models like the convexity 
of the region where the communication can happen [21]. Indeed, it seems to us that most of the arguments 
that we use are related to this property. 

The communication graph is given by the structure (V, E) where V is the set of nodes and E, the set of 
edges, i.e. the set of couples of nodes that can communicate together directly. We will use virtual edges. A 
virtual edge is an edge between two nodes x and y such that (xy) ^ E, but with a path from x to y of edges 
of E. 

Finally, we do not consider the impact of interferences or collisions during wireless communication. 



2.2 Coordinate system 

We use the virtual raw anchor coordinate system [4] with three anchors A\, A 2 , A3. It means that each node 
knows its distances to the three anchors, distances which form the node coordinates. I.e. the coordinates of 
node x is the vector (d(x, Ai),d(x, A 2 ), d(x, A3)). 

Definition 1. The coordinates of a node x is a vector 

(x 1 ,x 2 ,x 3 ) = (d(x,A 1 ),d(x,A 2 ),d(x,A 3 )) 

Throughout the paper, we suppose that all nodes lay inside the triangle defined by the three anchors on 
a 2D-plane, this area is denoted A. We use two different distances to define the coordinate system. 

In Section 3, we use the Euclidean distance for the distance function d. Given two points x and y, we 
note \xy\ the Euclidean distance from x to y, throughout the paper. 

In Section 4, we extend the results using for the distance d(x, Ax), d(x, A 2 ),d(x,A 3 ), the heights of the 
triangles A 2 xA 3 , A\xA 3 a.nd, AxxA 2 respectively. We will note this distance d (x, A t ) or d\. (x) for 1 < i < 3. 

We further suppose that AiA 2 A 3 is equilateral and that all nodes know the distances between the anchors: 

\A,A 2 l {AMlA^l 

3 Distributed graph planarization 

In this section, given an Unit Disk Graph G, we build a planar subgraph G. We further extend it to Gi by 
changing some of its edges by virtual edges, where a virtual edge represents a path of G. 

Recall that in [8] it is proven that if we consider a graph G — (V, E) and that we have three total order 
relations, <i, < 2l < 3 , on the set of nodes and 

— the intersection of the three order relations is empty, 

— for each edge (x,y) G E and for each vertex z ^ {x,y} there is at least one order <j such that x <j z 
and y <j z. 

then the graph admits a planar embeddings. In this paper, we adapt this result to UDG. It leads to a simple 
and localized distributed algorithm to planarizc a communication graph of a wireless network and to a simple 
description of the communication graphs that accept an efficient routing algorithm. 

Our aim is to define three suitable order relations by using virtual raw anchor coordinate system. The 
order relations have to satisfy some properties to ensure that the resulting communication graph admits a 
planar embedding. In this paper, we show how to locally compute the planar embedding by using only the 
distances to the anchors. We assume that there are no pair of nodes x, y such that for a given k G {1, 2, 3}, 
d(x,Ak) = d(y,Ak). It does not restrict the generality since the (Lebesgue) measure of these positions is 
zero. 

Given the nodes' coordinates, we define three total order relations, <i, <2> <3 on the set of nodes V in 
the following way: 

Definition 2. For k G {1,2,3} 7 nodes x and y with coordinates [x\ 1 x 2 ,x 3 ) and (3/1, y 2 , y 3 ) satisfy the 
relation x <& y, if and only if Xk < yk- 

Lemma 1. Given that the three anchors are not aligned, we consider the set of nodes that are inside A, 
see Figure 1. If we consider the restriction of the order relations <k on A x A denoted <k \a then their 
intersection is empty. 
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Proof. To prove that the intersection is empty is equivalent to prove that given any point x that belongs 
to the convex hull of the three anchors the triangular area A is covered by the three circles centered on the 
anchors and passing through x. Indeed, if the intersection is not empty there is a point y € A that belongs 
outside of the three circles (and reciprocally) , i.e. x <f. y, for k £ {1,2,3}. 

Because the area A is the union of the three triangles AxxA 3 , AixA 2 and A 2 xA 3 , see Figure 2(b), it is 
sufficient to show that the three triangles are covered by the circles. We consider AixA 3 particularly and 
the proof extend to the others triangles. We decompose the triangle A\xA-j, into two sub-triangles A\xxi 
and A^xxf, where xi is such that the line xxi crosses the line A1A3 perpendicularly. Because the length of 
the segment A\x is larger than the length of the segment A\xi the sub-triangle A\xxi is covered by the 
circle centered in A\ and passing through x. The same argument apply to the sub-triangle A^xxi and this 
concludes the proof. 




Fig. 1. If we do not restrict ourselves to the region A, (\=i <fc7^ 0- Dfc=i i s represented in gray 



Remark 1. Notice that if we do not assume that the nodes belong to the area A then, the intersection (1) 
may not be empty. Indeed there are point y whose the distances to the three anchors are larger than the 
distances of x to the three anchors, i.e. x < k y, Vfc = 1, 2, 3, see Figure 1. 



Definition 3. We define the three binary relations <i, <2, <3 by^x, y £ V, k = 1, 2,3, x< k y - 
k. 



x <k y and y <j x for j 



From Lemma 1, we deduce that the graph G~ c ~ a ~ r induced by these three total orders is planar. 

' ° ^ <1,<2,<3 J 

However, as we mentioned in the introduction, there are some major issues: 1) G~ c ~ v ~ T may not be a 
subgraph of an UDG, and 2) G~ ch ~ vd S r can not be computed locally. 

We denote mink the minimal z with respect to the order relation < k . The next lemma, gives a local 
condition to ensure planarity. 

Lemma 2. Given an UDG G = (V,E), if\/{x,y) £ E and Vz £ V \ {x,y} with max{|a;z|, \yz\} < \fhrj2, 
there exists k £ {1, 2, 3} such that x <k z and y <f c z then the graph is planar. 

Proof. The condition that max{|a;z|, \yz} < \fhrj2 is particular to UDG. Indeed, for an edge (u,v) £ E with 
either u or v at distance larger than \fhr jl from x and y, the two edges (xy) and (uv) (whose lengths are 
bounded by r) cannot intersect. This condition limits the set of nodes that potentially can be linked to an 
edge intersecting (xy) and ensures that the verification can be done locally. 

We now consider two edges (xy) and (uv), with both u and v at distance at most \/5r/2 from x and y. 
By assumption, there exist ki,k2,k3,k^ £ {1, 2, 3} such that 

u,v < kl x, u,v <k 2 y, x,y < ks u, x,y < ki v. 



It is clear that k\ ^ k%, k^ and &2 7^ k%, k^ and we can assume that k\ = ki and then u, v < / t 1 Tam(x, y). 
Indeed, if k^ ^ k^ we have k\ — ki because ki — 1,2,3. If k% — k^ we apply the same argument to u, v 
instead of x, y. 

We conclude that (uv) do not cross (xy) because each point of (uv) are <^ 1 smaller than x and y, see 
Figure 2(a). 
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Fig. 2. 



Lemma 3. Given a graph G = (V, _E) and three anchors Ai, A2, A3. We define the subgraph G — I V, E\ 

of G by Vx,y G V,(xy) G E <==> (xy) G E and 3k G {1,2,3} such that y = mirik{z \ x <k z} or 
x = mink{z | y <k z}. If the lengths of all the edges in the resulting graph G — I V,E) are smaller than 
2r/y5 ~ 0.8944 then the graph G is planar. 

Proof. Let (xy) be an edge of G. Since by hypothesis, all edges of G are of length at most 2r/v5, it can 
be seen as a subgraph of a Disk Graph with radius 2r/v5- To apply Lemma 2, it is sufficient to check that 
Vz ^ {x 7 y} with max(|aiz|, \yz\ < y/E/2 * 2r/v5 = r, there exists k G {1,2,3} such that x, y <k z. Only, 
if such a z exist, in our construction of G, we would not have the edge (xy), but instead, we would have an 
edge (xz), a contradiction. 

Remark 2. The selection procedure of the edges naturally induces an orientation. Hence the obtained graphs 
may be seen as digraphs. We will use this remark in Section 4. 

G is not a subgraph of G~ c ~ y ~ , however, Lemma 3 provides a sufficient condition ensuring that the 

subgraph G of the communication graph G is a planar graph, see Figure 3. The advantage of G is that 
it needs each node to know only its neighbors. It means that it is sufficient that each node broadcasts its 
identifier and its VRAC coordinates, so that its neighbors know them. It induces a communication complexity 
of 0(log(n)) bits. When the density is high enough, the proposed condition is sufficient. However, when the 
density is low, we observe that there are situations in which considering only the edges of length at most ^& 

leads to disconnect the graph. To avoid this, a solution is to reconstruct G~ c ~ y ~ r , but two questions arise: 
can we still do it locally ? If it is not a subgraph of the communication graph, how can we detect the missing 
edges ? For this, the solution that we propose is to introduce virtual links. If we are in the situation where 
x <k y and x <k z, z <f. y but z is out of range of communication of x, from Schnyder's theory, x would 
rather be connected to z than y. Only, due to the limited range of x, it does not occur. Then, the edge (xy) 
can potentially cross an edge from z if \zy\ < r. To avoid this, we replace (xy) by a virtual edge (xz). For 



this, the node y that knows its neighborhood informs x and a virtual edge between x and z (through y) 
replaces the edges (xy). In turn, z also checks if it is in the same situation as y. Ultimately, we would like 
that the computed graph Gi — (V, El) is a subgraph of the graph induced by the three total orders. 
The algorithm to compute the Gl goes as follow: 

— (As in Lemma 3) Each node x knows its neighboring nodes and compute the nodes yk, k = 1, 2, 3 such 
that yk = min{z | x <k z}. 

— (Virtual edge) Each node x checks with its neighboring nodes yp., k = 1,2, 3 that there does not exist 
a node yy in its second neighborhood such that j/&/ <k yk, % <k y 1 (d(yi,y) < r) and yy is out of the 
communication range of x. 

• If such a node does not exist the edge (xy) becomes active. 

• If such a node exists, yk check with yy that there does not exist a similar node that is out the range 
of communication of yi. This operation is repeated recursively until no node satisfying this property 
is found and a virtual edge is created between x and the last node found. The original edge (xyk) is 
removed. 

Remark 3. In the next section, we will see that, using the modified VRAC coordinates, a virtual edge 
represents a path of length 2. In this case, the recurrence is useless. However, using the original VRAC 
coordinates, a virtual edge can represent a longer path. 

In the middle of Figure 3 we plot a communication graph resulting from the selection described in Lemma 
3 without restricting the lengths of the edges to be smaller than -%. We observe that two edges cross. In 
the right of the Figure 3 the virtual links mechanism is used. We observe that the crossing is removed and 
the graph is planar. By comparing with the left side of the figure, we observe that the connectivity of the 
graph is better with the virtual links. However, when the node density is high enough, the selection of edges 
of length less than — % is sufficient. Furthermore, when the density increases, the number of virtual edges 
tends to zero (cf Section 6. 




Fig. 3. On the left the planar graph obtained by considering only edges of length at most -y=. On the middle the 
graph obtained by considering all the edges, we observe two edges crossing. On the right the planar graph obtained 
with virtual links (in green), the crossing is removed. 



Gi may not be a subgraph of G, which means that some of its edges may have length greater than r. For 
this, we can not use Lemma 2. However, if we can prove that it is a subgraph of G , we would obtain from 
[8], that it is planar. In the next section, we slightly change the VRAC coordinate system in order to prove 
that it is the case. It also allows us to give guarantees on the stretch of Gi. 

4 Properties of the planar embedding 

In this section we discuss a simple extension of the VRAC coordinate system. Using this new coordinate 
system, we prove that in G/, the distance between two nodes is at most twice the distance in the original 
graph G. Furthermore, we show that a virtual edge e corresponds to a path of at most two edges in G, and 
that the length of such an edge is upper bounded by 2r/V3. 
In this section, we make the following hypothesis: 

— There are three anchors A\, A2, A3, the nodes belong to the convex hull A of the anchors and they know 
their distances to all three anchors. 

— A1A2A3 is equilateral. 

— the nodes know the distances between the anchors (|j4i A2I, | j4i^4s| ,1^42^-3 1)- 

With respect to the first part of the paper, the two last hypothesis are new. By using the distances between 
the anchors, each node x can compute the heights of the triangles A2XA3, AixA^ and, AxxA^. As A1A2A3 
is equilateral, it is equivalent to compute the surface of these triangles or the heights and it is then easy to 
see that their sum is constant. We denote these values (xi,X2,xs), see Figure 4(a). One advantage of this 
coordinate system with respect to just using the distances to the anchors is that the sum of the three triangle 
areas is constant, so we can normalize the coordinates such that x\ + X2 + £3 = 1- A reason for doing this is 
that because of the measurement errors on the physical location of the nodes, it is likely that the distances 
to the anchors do not correspond to coordinates inside a same plane. With the normalization, we project 
the coordinate on a same plane. 

4.1 Adapting results of Section 3 and further 

Results of Section 3 Using the coordinates defined above, we define the order relations <i,<2,<3 and 
<ii <2j <3 the same way we did in Section 3. In Section 3, given a node x, the nodes satisfying y>kX were 
outside the circle centered at Ak of radius |xA^|. With the new definition of the distance function d (c.f. 
Section 2.2), the nodes y satisfying y>kX are contained on the half plane containing Ak defined by the line 
parallel to (Ak mod 3+iAk mod 3+1) going through x, as illustrated in Figure 4(b). Using this observation, it 
is easy to see that the intersection of the three order relations is empty, so Lemma 1 is still valid in this 
coordinate system. Similarly, Lemma 2 remains also true and we can adapt the proofs of Lemma 3. 

To summarize, all the results we have proved previously are valid with the new coordinate 
system. 

Connectivity results and stretch 

Definition (Figure 4(d)): Given a node x, we call the greedy regions of x the three regions Af = {z \ 
x <i z}, for i e {1,2,3}. 

Definition (Figure 4(d)): Given a node x, for i G {1, 2, 3}, we denote the region between the two regions 
^and^ mod3+1 byAf. 

Remark 4- A node x has at most one outgoing edge towards a node in each of its greedy regions. It has no 
outgoing edge towards node not in its greedy regions, however, it may have an ingoing edge from any node. 




(a) The new coordinate system: 
x — (xi,X2,xs) where xi, xi and 
xz are respectively the heights 
of the triangles A2XA3, A1XA3 
a,nd,AixA2- 




A, 




(b) y with x< 3 y. 




(c) New proof of Lemma 2 

Fig. 4. Adapting proofs to the new settings 



(d) Definition of the greedy re- 
gions 



Lemma 4. Given an edge (x, y) £ E, there is a path P from x to y in Gl. The path is contained in either 
{z £ Af\z <i y} or {z £ A\\z <, x} for some i £ {1,2,3}, and it verifies: X) e pp l e l — 2|a?y | . 

Proof. Given an edge (x, y) £ E, without loss of generality, we can suppose y £ Af. (x, y) £ E. By hypothesis, 
\xy\ < r. 

We prove the lemma by induction on the length of xy 2 . By our definition and Lemma 3, (x, y) £ El iff 
y = jnin\{z \ x <\ z}. 

If y = min\{z \ x <\ z}, then (x,y) £ El and there is a (direct) path between x and y. 

If not, there is x' with x'<i y and x' — min\{z \ x <\ z}. Notice that x' £ {z £ A x t \z <j y}. 

We now prove that \x'y\ < \xy\. 

Using the notation of Figure 5(a), we obtain that \x'y\ is maximum for x' = x, x' = c\ or x' — cy.. We 
hence have to prove \c^y\ < \xy\. The case \c\y\ < \xy\ is symmetric. 

In the purple triangle of Figure 5(a), we have sm ^%_ p) = si ^ /3) . So gf = ""ff^^ ■ As /3 > tt/3, 

we have ^ < 1. So \c 2 y\ < \xy\. 

Recall that we supposed that no two nodes in the network have the same coordinates compared to a given 
anchor. It means that x' can neither be C2 nor ci 3 . We hence have \x'y\ < max(\xy\, \c\y\, \c%y\) < \xy\ < r. 

So we can apply the induction on the edge {x'y) which is in E and strictly shorter than (xy) and this 
proves that there is a path between x and y through x' . 



2 The length of an edge takes value in 1R, however there are a finite number of edges (upper bounded by n 2 ), hence 

our induction will terminate. 
3 



Because d(x, A2) — d(ci, A2) and d(x, A3) = d{c2,Az) 



A, 




(a) 




Fig. 5. \x'y\ < \xy\ and E eeP l e < 2 \ x v\\ 



Let us now look at the stretch factor. We have x <i x' , y >i x' and there is a path y 1 — x' ,y 2 , . . . ,y l = y 
in the graph G that satisfies y 1 <i y 2 <i . . . <i y l — y by construction of the virtual edge. Because the 
coordinate with respect to A\ increases(monotonically) along the path we have that \x\ 



\y x — y[\ = \xi ~ 2/1 1 (the subscript indicates the coordinate with respect to A\). 



y\\ 



\vi 



y\\ 



If x' £ A\ for i £ {2, 3}, we have y <t x' and y >i x' . By induction, if such a i exists, the rest of the path 
will be in A\ and the ith-coordinate decreases (monotonically) along the path. If not, x' £ A\ and y >i y' , 
y <2 x' and y < 3 x' and both coordinates with respect to A 2 and A^ decreases along the path. This proves 
that in all cases, we have that there exists i £ {2, 3} such that \xi — y}\ + \y} — yf\ + .. ■ + \yi~ 1 — y\\ = \%i— y%\- 

In summary, along the path P, for z £ P, z\ increases from x\ to y± and there is an i £ {2, 3} such that 
Zi decreases from Xi to yi. The distance covered in A\ to go from a node with i th coordinate xi to a node 
with i th coordinate yi is upper-bounded by —/=\Xi — yi\, see Lemma 5, since we suppose that the anchors 
form an equilateral triangle. 

From this we deduce that the length of the path P is upper bounded by -t= \xi — yi\ + -j^\%i ~ yi\, 
this longest path is obtained by moving along the path where the ith coordinate is constant and the first 
coordinate goes from x\ to y\ and then along the path where the ith coordinate goes from Xi to yi and the 
first one is constant. 

We now express \xy\ in terms of \x\ — yi| and \xi — y%\. By the configuration of the different triangles, c.f. 

]2 +{ J rM-y l \- J r,W-yi\) 2 

-vd+^ki- 



Figure 5(b), we have \xy\ 
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Corollary 1. Gl is a subgraph of G which is equal to the half-O^ graph. 

Proof. The definition of the half-6>6 graph gives the same graph as G when using the three total orders 
using the modified VRAC coordinates. From the previous lemma, we get that if there is an edge {xy) £ El, 
with y in a greedy region A^, then y is minimum according to <&. Hence, (xy) is an edge of G s . However, 
there are examples in which Gl 7^ G 

From this corollary, we immediately obtain that Gl is planar. When Gl = G s , we also deduce that it is 
geometric 2-spanner. However, in the general case, we can not deduce any information on its stretch, which 
is the object of the next theorem, which is implied by Lemma 4. 

Theorem 1. Given a connected graph G, the graph Gl is planar, and for any two nodes x and y, if there is 
a path of length (. from x to y in G, there is one of length at most 2£ in Gl. 
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However, notice that the previous theorem applies to Gi which may contains virtual edges that are not 
edges of the Unit Disk Graph G, instead, a virtual edge represents a path in G. The next lemma says that 
such a path has length two. Notice further that the virtual edges are edges of the Unit Hexagonal Graph 
(c.f. [16]). Finally, simulations in Section 6 show that the virtual edges are rare and disappear as the node 
density increases. 

Lemma 5. A virtual link (xy) G El has length at most 2r/v3 when the anchors form an equilateral triangle. 

Proof. Without loss of generality, we can suppose that the link is oriented from x to y and with respect to 
minimizing coordinate y\. The virtual links represents a path x, z, . . . , y where z is inside the communication 
range of x, so d(x, z) < r. We know by construction that considering the first coordinates, we have 2/1 < Z\. 
It means that y is in the triangle T\ delimited by the three lines \u\u\ — zi}, \y\v2 = x-i\ and {u>|w3 = X3} 
as depicted in Figure 6. The furthest points of this triangle are the two summits other than x. This triangle 
is equilateral and the edges are of length 2r/v / 3, so \xy\ < 2/v / 3r. 

From this lemma, we obtain that Algorithm 1.1 constructs correctly Gi. 

Input: A Unit Disk Graph G. 
Output: Gi. 
for all x G V do 

for k G 1, 2, 3, y k - min k {y G A%, \y k x\ < r} do 
x broadcast "activate (xy k )" 
y' k = min k {y G A%, \y k y\ < r} 
if yk ■£ y' k then 

y k broadcast "disable (xy k ) and activate (xy' k )" 
end if 
end for 
end for 

Algorithm 1.1. Distributed construction of Gi. 



Corollary 2. Given a connected graph G, if the density is high enough, Gi is a planar 2-spanner of G. 

Proof. When the density increases, the number of virtual edges decreases, and if the density is high enough, 
with high probability, there remains none of them. In this case, Gi is a subgraph of G, and the result follows 
from Theorem 1. 

To summarize, Gi is a planar graph such that 1) the length of a shortest path in this graph is at most 
twice the length of a shortest path in the communication graph, and 2) virtual edges correspond to a path of 
length two in the communication graph. Furthermore, to construct Gi (Algorithm 1.1), the communication 
complexity in terms of bits at a node x is as follow: each node broadcast once its Id and coordinates, then, 
each node broadcast the Id of its three neighbors minimizing the orders, and finally, each of this neighbors 
may send to x an other node Id which is the extremity of a virtual edges starting from x. Hence, each node 
induces an exchange of at most 6 nodes (excluding the neighborhood discovery), which gives a total of at 
most 6n nodes' ids that are broadcast. The computational complexity is 0(nA), with A the maximum degree 
of the communication graph: each node x computes which of its neighbors minimizes each order, plus, x 
requires each of the selected neighbors 2/1,2/2, 2/3 to verify if there should be a virtual edge, which also consists 
in computing a minimum of a set of at most A elements. 

Now, if instead of using VRAC coordinates, we use the Euclidean coordinates, each node y can computes 
on its own if it is minimum for one of the three orders for one of its neighbors. It means that we can avoid 
the statement "broadcast "activate (xj/fc)"" in Algorithm 1.1. The modified version of Algorithm 1.1 needs 
the broadcasts of at most in nodes identifiers that can all be performed in a single round of communication. 
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In case there arc no virtual edges, the resulting graph is a 2-spanner of the Unit Disk Graph, and in this case 
no messages are exchanged. It answers the open question 22 of [16] under the hypothesis that the density id 
high enough. Recall that Gl is always a subgraph of a Unit Hexagonal Graph (c.f. [16]), and hence a planar 
spanner for these graphs. 

5 A local routing algorithm 

We now propose a local routing algorithm. This algorithm has two modes depending on if the destination is 
in a greedy region of the sender or not. 

Lemma 6. Recall that we suppose the triangle AiA 2 A% equilateral. Let x be a node with a message for a 
destination y. If y belongs to a greedy region of x, and that x has an out-neighbor in this greedy region, the 
algorithm proceeds as follow: 

— (Data delivery) \xy\ < r in which case x transmits the data directly to y. 

— (Greedy routing) \xy\ > r and x transmits to its neighbor x' that belongs to the same greedy region as 
y. We have \xy\ > \x'y\. 

Proof. If x transmits directly to y there is nothing to prove. For the other case, let x' be the node that 
receives the message. We note the vector xy — ae lai and xx' = fee M2 (\xy\ — a and a > b). The vector 
x'y = ac tai — fee 1 " 2 , \x'y\ = a 2 + b 2 — 2abcos(ai — a 2 ) and, < ai, a 2 < a because x' and y belong to the 
same greedy region. Then, —a < ct\ — a 2 < a and cos(ai — a 2 ) > 1/2. This is sufficient to conclude that 
\x'y\ < \xy\ =a. 

Lemma 7. We assume that the anchors form an equilateral triangle. Let x be a node that wants to send 
a message to a destination y using the graph Gl or direct transmission. We assume that y does not belong 
to a greedy region of x. Without loss of generality, y is in A 2 . We further assume that x has three out- 
neighbors and so do any node in the equilateral triangle T with base the segment parallel to (^42^.3)7 centered 
at x, of length 4/y3r and with the other summit in A 2 (c.f. Figure 6). Under those hypothesis, we have 
two paths (without consideration on the orientation of the edges) P\ — x,uq, P® , Mi, ..., Uk-i, Pi , u k,z and 
P 2 = x, Vq, P 2 , v\, ..., vi-i, P 2 ~ ,Vi,z, without loss of generality uq >i vq (when P 2 7^ x,z, as if P 2 — x,z, 
there are no vq), and we have: 

— k- 1 < I < k. 

— the P? , for i G {1,2} and j < k are monotone paths with respect to > 1; potentially of length 0. 

— VuePi\{z}, us A%. 

— VveP 2 \{z}, v€A%. 

— VO < i < I, there is an oriented edge from m to m and, \/u £ {uj, PI}, u >i Vi- 

— VO < i < I — 1, there is an oriented edge from Vi to Uj+i and, Vw G {vi, P 2 }, V >i V»+i. 

— z is in A 2 . 

Given these two paths, either a node from {x, u , ...Uk, v 0l ...Vi, z} has y within its communication range, 
or \zy\ < \xy\. 

Proof. By hypothesis, x has a neighbor (in Gl) uq in A 2 and a neighbor vq in j4|, without loss of generality 
uq >i vq. By Lemma 5, which bounds the length of an edge, uq is in T, so uq has three out-neighbors. In 
particular it has an out-neighbor v! in A^° . u' can not be in A^° ("1 A 2 as otherwise we would have a virtual 
link (xv!) instead of the link (xuq). 

If v! is in A 2 , if we can prove that (xv!) G El setting z — v! ', the theorem would be verified with 
Pi = x, uq, z and P 2 = x, z. But, we may not have an edge between x and v! , however we prove that u' can 
be replaced by a z" satisfying all the desired properties, v! is in T, so it has an out-neighbor u" in A\ . If 
u" 7^ x, u" has an out-neighbor in the greedy region oriented towards A 2 that in turn has a neighbor in its 
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Fig. 6. On the left side, the three greedy regions associated with x and the edges to the three nodes u — miri2{z 
x<2 z}, v — mins{z \ x<3 z} and, w = min\{z \ x<i z}. On the right, a path leading to z. 



greedy region oriented towards A2, ... By planarity, since we cannot cross the edge uqz this path has to go 
through vq. So we have a path from u" to uq. Similarly, u" has an out-neighbor u" in A\ which is either 
x or has a path leading to uq- Ultimately, we will have a node z' with out-neighbor x and a path to Uq. 
Using the same argument and considering neighbors in A\ , we have that z' has u has neighbor in A\ , or 
it has w whose out-neighbor in A™ is x, and that has a path leading to uq in A^. Going along this path, we 
reach z" whose out-neighbor in A\ is x, and out-neighbor in A\ is u§. Hence the theorem is verified for 
P\ = x, uq, z" and P2 = x, z" 

We now suppose that v! is in A| and we want to prove u' — vq. By contradiction suppose that u' =^ v . 
Recall that we have u' <i uq, by construction of Gi as well as u' < 3i»o- By planarity, we have vq >i u' . 

The node u' has a neighbor in its greedy region A\ . This neighbor must be ^o- Indeed, if not, using 
v! < 3vq and using the same argument as before, this neighbor it" has a neighbor in the greedy region A^ 
that in turn has a neighbor in its greedy region oriented towards A2, ... But, since we cannot cross the edge 
uqu' by planarity, there is an infinity of such nodes which is impossible. So u' has vq has neighbor in its 
greedy region A l { . But then vq has a neighbor in its greedy region A"^ . Looking for a sequence of neighbors 
in the greedy region oriented towards A2, will ultimately lead to cross edge (uqu') which contradict the 
planarity of Gi. 

So v! = Vq 



Since vq is in T, it has a neighbor v' in its greedy region oriented towards A2. v' cannot be in A^° , the 
same region as v because if it was true, by construction of G/, there would be an edge between m and v' 
instead of the edge (uqVo)- 



First suppose that v' € Af . We would like to set z — v' , but, as for the very first case considered in 
the proof, it may not be appropriate. Indeed, there may not be an edge (zvq). So we proceed as follow: 
we consider the out-neighbor u" of v' in A\ . If u" = vq, then we are done by setting z — v' . Else, we can 
prove that if has Uq as out-neighbor in A^ . We proceed recursively until we reach i>"' which has Uq as 
out- neighbor in Af and vq as out-neighbor in A\ (Notice that uq can not have v" ' as out-neighbor in A^° 
as it already has vq). We set z = v"' and the theorem is proved. 
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We now study the case v' G A\ ■ Let call v" the out-neighbor of v' that is in the greedy region oriented 
towards A\. If t>" = t>o, all is fine. Else, we must have v" >o u (otherwise the neighbors directed towards 
A3 will cross the edge vqv'). Considering the following out-neighbors in the greedy region oriented towards 
A\, we will reach a node v"' whose out-neighbor in A\ is Uo, or which is the out-neighbor of x in A%. We 
call this path from v" to v"', Pf and we set v' = v\. 



We continue similarly by considering the out-neighbor of v± in A\ 



To see that this process terminates, notice that the first coordinate of the vertices on the paths we 
construct decrease as we get closer to the line (^42^3), plus the distance between a node in A\ and a node 
in .A3 is lower bounded by a bound which increases when we get closer to the line (A2A3), plus the length 
of an edge is upper bounded by Lemma 5. Hence the process has to finish on a vertex z in A%. 

We now prove the second part of the Lemma. The maximum length of a virtual edge is 2/v / 3r, hence we 
have that the polygon formed by x, u , ...Uk, z,vi, ...Vq is composed of triangles with two side of length at most 
2/\/3- So it is covered by the union of the disk of radius r centered on {x, u a , ...Uk, i>o> ■■■ v h z}. If V is inside 
the polygon formed by x, uq, ...Uk, z, vi, ...vo, it is within the communication radius of one of these vertices. 
If not, y is bellow the polygon x, Uq, ...Uk, z, vi, ...vo and out of range from both vi and Uk in particular, and 
then \zy\ < \xy\ . 

Theorem 2 (Zig-Zag : an extended greedy routing). 

We assume that the anchors form an equilateral triangle. Let x be a node that needs to transmit a message 
to a destination node y. If any node at distance less than 4/v3r of x has three out-going neighbors, then the 
following strategy delivers the data either to y or to a node z closer to y than x. 

— If y is in the communication range of x, x sends the message to y. 

— If y is in a greedy region of x, x sends the message to its out-neighbor which is in the same greedy region. 

— Otherwise, use the restricted greedy routing process starting at x. Wlog y G A^. 

• x sends the message to its out-neighbor in A^ U A% which has the highest first coordinate. 

• A node u G A% sends the message to y if possible or to its out-neighbor v verifying u >i v and 
v >3 x. If v >3 x and v >2 x (i.e. u G A^), end the restricted greedy routing process. 

• A node v G A| sends the message to y if possible or to its out-neighbor u verifying v >i u and 
u > 2 x. Ifu> 3 x and u > 2 x (i.e. u G A^), end the restricted greedy routing process. 

Proof. By applying Lemmas 6 and 7 we see that the routing strategy leads to y or to a node that is closer 
to y than x (z in Lemma 7). Indeed, Lemma 7 ensures that a restricted greedy routing process starting at 
x follows the path uo,Vo,Ui, ..., z using the same notations. 

6 Simulations 

We implemented both the planarization algorithm and Zig-Zag, the routing algorithm, and we present bellow 
the results of the simulations. We considered a network composed of 300 sensors spread in a square area 
[0; 1] x [0; 1] with three anchors at position (0.5, 3.5), ( — y= +0.5, —1.5) and (-7= +0.5, —1.5). We considered a 
communication radius r for the sensors which ranges from 0.11 to 0.225. For each value of the communication 
radius, we performed the average over 1000 networks and successfully routed messages. Notice that the value 
are plotted with respect to the average degree of the nodes in the UDG, as in G7, the average degree is upper 
bounded by six since it is planar. 
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We first plot the number of virtual edges in Figure 7(a). The simulations indicate that the number of 
virtual edges tends to when the node density increase. Indeed, the average number of virtual edges decreases 
from 1.6 when the network is sparse, to 0.03 when it is dense. 

Then, in Figure 7(b), we plot the average stretch of the path computed in Gl by Zig-Zag, where the 
stretch of a path is the length of the computed path divided by the euclidean distance between the source 
and the destination. We observe that we have a stretch which is between 1.3 and 1.4 which is better than 
the theoretical stretch factor of 2, and we compare it to the stretch of the greedy algorithm (in Gl), which 
is slightly better and of approximately 1.3. This apparent gain of efficiency is to be mitigated by the results 
shown on Figure 7(c) which indicates the success rate of the greedy algorithm and of Zig-Zag. The success 
rate of Zig-Zag tends to 100% as the density increases, which backen the theoretical results of Theorem 
2. Indeed, when the density increases, the hypothesis of Theorem 2 arc verified, whereas they may not be 
verified at some nodes when it is low, thus reducing the success rate. The success rate of the greedy algorithm 
also increases with the node density, however, it remains bellow 80%. 
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